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Abstract 

We investigate the G/G gauged Wess-Zumino-Witten model on a Riemann sur- 
face from the point of view of the algebraic Bethe Ansatz for the phase model. 
After localization procedure is applied to the G/G gauged Wess-Zumino-Witten 
model, the diagonal components for group elements satisfy Bethe Ansatz equations 
for the phase model. We show that the partition function of the G/G gauged Wess- 
Zumino-Witten model is identified as the summation of norms with respect to all 
the eigenstates of the Hamiltonian with the fixed number of particles in the phase 
model. We also consider relations between the Chern-Simons theory on 5 1 x S/ 8 
and the phase model. 



1 Introduction 



It is known that there exist connections between topological gauge theories and the Bethe 
Ansatz for integrable systems (Gauge/Bethe correspondence). In [I], the equivariant lo- 
calization for the BF theory coupled to a one-form valued adjoint Higgs field (Yang-Mills- 
Higgs system) is studied and the localization configurations lead to Bethe Ansatz equa- 
tions for the non-linear Schrodinger model. Further, the partition function of the Yang- 
Mills-Higgs system is related to norms of wave functions for the non-linear Schrodinger 
model [2]. So we expect that the partition functions of other topological gauge theories 
are also related to norms of wavefunctions for the corresponding integrable systems. 

Wess-Zumino-Witten (WZW) model of a two-dimensional conformal field theory has 
rich structures and various applications in mathematics and also in physics. For example, 
the Hilbert space of the Chern-Simons (CS) theory with a gauge group G on R x E/, is 
equivalent to the space of the conformal block for the G/G gauged WZW model on a 
Riemann surface E^. The partition function of the CS theory on a three manifold can 
be obtained by sewing the boundary Riemann surfaces which is implemented by an inner 
product of states on E. One can also calculate Wilson loop expectation values which give 
knot invariants in terms of fusion coefficients and modular matrices [3]. 

In the WZW model, one can construct the G/H gauged WZW model by gauging an 
anomaly free subgroup H of the global symmetry group G. The G/H gauged WZW 
model is an explicit lagrangian realization of the coset construction in the CFT. When 
H = G, the G/G gauged WZW model becomes a topological field theory [3], [5]. One 
expects that there exists a method for calculating the partition function and correlation 
functions without relying on the CFT techniques nor the representation theory of the 
affine Lie algebra. Actually, these were derived by a field theoretic approach in [6], [7]. 
In this approach, it is important for the G/G gauged WZW model to possess a certain 
BRST-type symmetry whose square generates a G-gauge transformation. This symmetry 
makes it possible to work out the path integrals with insertions of BRST closed operators 
via equivariant localization procedure. In higher rank of the gauge group, the localization 
configurations for the diagonal components of G-elements are complicated. However, the 
final expression for the partition function is simply expressed by modular matrices and 
an integrable structure might inhibit behind the localization configurations. 

In this paper, we will show that the integrable system corresponding to the U(N)/U(N) 
gauged WZW model is the phase model. The phase model is an integrable field theory 
on one-dimensional lattice. We can apply the algebraic Bethe Ansatz. For example, see 
[8], [9]. It is known that the phase model appears in the SU(N) WZW model. Recently, 
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Korff and Stroppel established the su(N)k Verlinde algebra [10] in terms of the algebraic 
Bethe Ansatz for the phase model and derived an efficient recursion relation for calculat- 
ing fusion coefficients [11] . See also a short review [12]. We will consider relations between 
the Gauge/Bethe correspondence and [TT]. We also point out that the partition function 
of the CS theory on S 1 x E^ is related to norms of Hamiltonian eigenstates for the phase 
model. 

This paper is organized as follows. In section |2j we study the localization for the 
U(N)/U(N) gauged WZW model and find that the localization configurations satisfy 
particular equations. In section [31 we review the algebraic Bethe Ansatz for the phase 
model. In section HI we find that the Bethe Ansatz equations equal to the localization 
configurations of the U(N)/U(N) gauged WZW model under parameter identifications. 
We find relations between the partition function of the U(N)/U(N) gauged WZW model 
and the Bethe norms for the phase model. The section [5] is devoted to conclusion. 

2 G/G gauged Wess-Zumino-Witten model and lo- 
calization 

In this section, we calculate the partition function of the G/G gauged WZW model on a 
genus-/i Riemann surface E^ by applying localization method. This model is a topological 
field theory and its partition function counts the number of conformal blocks in the G 
WZW model with level k on E^. When we set g = e l ^/ k and expand them at the leading 
order of 1/k, the G/G gauged WZW model on T, h reduces to the BF theory with the 
gauge group G on [13]. This can be regarded as an another deformation of the BF 
theory different from the Yang-Mills-Higgs system. 

The partition function of the G/G gauged WZW model on E^ is defined by 

^GWZw(Sh) = J VgVAV\e- kS ^> A > x \ (2.1) 

The action of the G/G gauged WZW model is 

<Sgwzw(0, A A) = S wzw {g) ~ 7T I zYx{A z d- z gg~ l - A z g~ l d z g - g^A z gA z + A Z A 2 ) 
__L f d 2 zTr(\ z \ 2 ) (2.2) 

with 

Swzw(g) = -^ [ d 2 zTi{g- l d z g-g- l d z g)-iY{g) (2.3) 
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and the Wess-Zumino term T(g) is 

r(9) = 7^ [ d^Tr {g- l d l9 ■ g-%g ■ g- l d k g) . (2.4) 
12tt J b 

Here B is a three dimensional manifold with the boundary dB = £/,,. g(z,z) is the G- 
valued field, A = A z dz + A z dz is a two-dimensional gauge field and A = X z dz + X z dz is 
a one-form adjoint fermion. We also denote the holomorphic part of A as A^ 1 ' ^ = A z dz 
and anti-holomorphic part as A^ ' 1 ^ = A z dz and so on. Q is a scalar BRST charge defined 
as 

QA = A, QA (1,0) = (A*)™ - A (1,0) , QA (0 ' 1} = -(A* _1 ) (0,1) + A^\ Qg = 

(2.5) 

with A 9 = g~ 1 dg + gAg^ 1 . The partition function (12. ip is invariant under the BRST 
transformation. The square of the BRST transformation generates following gauge trans- 
formations 

C g A^ = (A 9 ) (1,0) - A^'°\ CgA^ = -(A^fV + A m t 

C a \W> = g-WWg - A (1,0) , £ 3 A (0,1) = -gX^g- 1 + A (0,1) , £ g5 = 0, (2.6) 

where Q 2 = C g . 

From now on, we set the gauge group G as U(N) for simplicity. We evaluate the 
partition function (I2.ip . The resulting expression is 

1 oo r N N 

ni,--- ,rt]v=— oo a=l a,6=l 

xexp|27rzEn Q ^(iV + A;)0 a - £ b + j J . (2.7) 

Let us briefly explain the derivation of (12. 7p . See detailed calculation [jBj. First of all, 
we take a diagonal gauge on g(z, z) G U(N) at (12. ip . that is, g(z, z) = e 2m ^ a = 1<f ' aHa where 
H a are the Cartan generators. Then the path integral with respect to g(z, z) becomes 

TV 

! Vgf(g) = —f Y[V4> a Bet crc (l - ad(e 2 ^))/(e 2 ^), (2.8) 

J ' J a=l 

where \W\ = N\ is the order of Weyl group for the gauge group U(N). Here Det Ci g is the 
Faddeev-Popov determinant for the diagonal gauge fixing. Next, the path integral with 
respect to off-diagonal components of the gauge field gives a functional determinant 

Det- 1 1 , 0)A(0 , 1) (l-ad(e 2 ^)). (2.9) 
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The each functional determinant of the ghost and the gauge field are infinite dimen- 
sional and divergent, but their differences can be explicitly evaluated by the Hirzebruch- 
Riemann-Roch theorem: 



Det C|g (l ad(e 27 " )) _ TT n _ 2n%(4> a -4> b )\i-h+n a -n b / 9 in s 

Det A(li o )jA( o, 1) (l-ad(e^)) " > ' ^ 

a^b 

where n a is an a-th diagonal [/(l)-charge of the background gauge field: 

n a = -L [ fW. (2.11) 

Combining these results, the partition function of the U(N)/U(N) gauged WZW model 
is 

N N 



I VgVAV\e~ kSGWZW = ^ I f[ Vcf> a VA {a) VX f[ (l - e 



xe W (ikJ2 I AAA )' ( 2 - 12 ) 



a^b 



a=l " ^ h 

Thus, the effective abelianized action obtained the BF-type action 

A? 



^ / <P a F a . (2.13) 



o=l 

Here we note the fermion bilinear term AAA. Since the effective action (I2.12p is 
not BRST invariant, we have to add appropriate counter terms to restore the BRST 
symmetry by requiring the effective action to satisfy descent equations. This leads to 
level-shift k — >■ k + iV for the coefficient of the fermion bilinear term. Of course, the 
such renormalization do not makes influences on the later calculation because the fermion 
bilinear term enters in the effective action freely. The renormalization effect becomes 
crucial when we couple the theory with additional matters [1] |14j . 

F b can be always decomposed to a harmonic part and an exterior derivative of a 
one-form da^ such that 

F h = F {h) +da b . (2.14) 

Integrating by part puts delta functional constraints on d<f) a , so the fields <f) a (z,z) 
reduce to constant fields. Therefore, we obtain (12.71) . 
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We rewrite (I2.7P by using the Poisson resummation formula: 

1 oo „ N N 

«w(s.) = ^ e /n^ n a - e 2 *-^) 1 -* x 

mi,- ,mjv=— oo a=l a, 6=1 



a^6 

N / A 

X 

a=l \ 6=1 



n ^ ( (iv + fc)0 a - e ^ + - ma ) • (2 - is) 



The partition function (I2.15P is invariant under the interchange k <-+ N because the 
U(N)/U(N) gauged WZW model on has a property of the level-rank duality [T5] . 
Therefore, we can rewrite (I2.15P as 

oo n k k 

zsslph) =4 e / n ^ n c 1 - e 2 *-^ ) ) i - ft x 

m ir"i m t =_ °° a= l a, 6=1 

a^6 

fc / fc , _ , \ 

II 5 (JV + A;)0 a - E <Pb + —z- ~ m a . (2.16) 



x 

a=l \ 6=1 



Integrating (I2.16P with respect to 4> a s, the partition function localizes to configurations 
which the constant fields cj) a satisfy constraints 



(iV + A:)0 a -E^ + ^--^a = O. (2.17) 

6=1 

Let us consider solutions of these equations. We immediately find that the solutions 



are 



' f Ja+h\J\\), (2.18) 



™ k + N \ N 

where J a = m a — m a G Z and || J\\ = Yjo=i for a = 1, ■ ■ ■ ,k. Note that range of 
the each field <fr a is < <ft a < 1 . The fields in this range only contribute to the partition 
function. Therefore, we count the number of solutions of these equations in the range 
< (j) a < 1. We immediately notice that when (fi a = fa for a ^ b the configurations do 
not contribute to the partition function. All the <p a are contained in the range, even if 
we interchange the all solutions <p a . So we can set 0i < 02 < • • • < <Pk and a factor A;! in 
the partition function cancels out. The number of piecewise independent solutions of the 
equations in the range of < <p a < 1 is 

(jfe + JV-1)! , 
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This number coincides with the number of primary fields of the SU(N)k WZW model. 
The each solution (" 12 . 1 8[) is in one-to-one correspondence with the primary fields of the 
SU(N)k WZW model or the highest weights of the integrable representation in the afline 
Lie algebra siifc(iV). When the set of this solutions is denoted by {Sol}, the partition 
function is 

fn'<-i(e 2 "*"-e 2 ™**)1 1 ~" 

</»!,-, fc £{Sol} I lla=l e I 

where a and (3 are a genus independent and dependent constant, respectively. 

Finally, we determine the normalization for the partition function of the G/G gauged 
WZW model on which is compatible with the number of the conformal blocks in the 
G WZW model on E/j. The partition function of the G/G gauged WZW model on 
also can be represented by the modular S-matrix for the character in the G WZW model 
as follows 

^gwzw(^) = y~l(^o^) 2 2h ■> (2-21) 

TZ 

where TZ denotes an integrable highest weight representation in the afline Lie algebra q 
corresponding to a primary field in U(N) WZW model and the summation runs through 
all the primary fields in the G WZW model. Therefore, we determine the normalization 
such that the partition function (I2.23f) matches with (I2.2ip in G = U(N). When h — 1, 
the partition function of the G/G gauged WZW model coincides with the number of 
the primary fields in G WZW model. The genus independent normalization factor a is 
(N + k)/N because {Sol} only runs through the primary fields in the su(N)k WZW model 
and the number of the primary fields in the ii(N)k WZW model is 

(k + N)\ 
N\ k\ ' 

The resulting partition function of the U(N)/U(N) GWZW model on is 

N + k j 1 Ui b =i(e 2 ^-e^] sl - h 

4wzw(^) = ^ r E i (k + m n k ew°-w° ' ' (2 ' 23) 

3 The phase model 

The phase model is a quantum integrable field theory on one-dimensional lattice and a 
strongly correlated boson system. We can calculate observables such that energy eigen- 
values, eigenvectors and correlation functions and so on in this model, by applying the 



(2.22) 
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algebraic Bethe Ansatz to this model. This model is considered by e.g. [8], [9], [TT]. In 
this section, we review the phase model and the algebraic Bethe Ansatz for this model. 

3.1 The phase model 

In this subsection, we define the phase model. First of all, we define the phase algebra. 
The phase algebra $ is an algebra such that operators {N, (p, obey 

[N,<p] = -<P, [N,<p*]=<p\ W f = l. (3.1) 

The operators ip and ip^ serve as an annihilation operator and a creation operator, re- 
spectively. Next, we define a Fock space T for the phase algebra given by the equations 
(ETIj) . A vacuum state |0) is defined as (p\0) = 0. The set {|m) := (^ f ) m |0) | m G Z> } 
forms the basis of the Fock space by successive actions of the creation operator p^ on the 
vacuum state |0). It also holds following relations, 

N\m) = m\m), ip^\m) = \m + 1), p\m) = \m — 1). (3-2) 

The Hamiltonian of the phase model on the one-dimensional lattice with the total site 
number N is given by 

1 N 

H = -- ^ (fM+i + <Pi<Pi+i) (3-3) 
i=i 

where we imposed a periodic boundary condition N + 1 = 1 and set the lattice spacing 
A = 1. The operators {p>i, <p\, iV* }*=!,... ,jv obey the iV-fold tensor product $ 0JV of the 
phase algebra (13. ip . $® is defined by 

ViVj = <Pj<Pi, <p\<Pj = NiNj = NjNi 

<Pi<Pi = 1 , <Pi<p] = <Pj<Pi if i ^ 3 

N l {i- v y l ) = Q = (i-^y l )N t . (3.4) 

where indices i,j label the sites of the lattice. Therefore, the Hamiltonian belongs to <&® N 
and acts on the iV-fold tensor products of the Fock space J 7 ®^. The basis of J 71 ^ consists 
of {\m 1 , ■ • • , m N ) := \mi) <g) • • • ® \m N ) \rrii G Z> }. 

To better understand this model, let us change the operators obeying the phase algebra 
to the operators {ai,a\,Ni} i= i t ...^ obeying the free boson algebra 

[Ni, dj] = Sijdj, [Ni, a]] = [a i ,a^] = 5 ij , Ni = a\ai, (3.5) 
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as follow 



Vl + Ni 



(3.6) 



where (1 + iVj) -1 / 2 is defined as formal power series. Substituting (I3.6P into the Hamil- 
tonian (I3.3p . we found that the Hamiltonian has infinite interaction terms in front of the 
hopping term. Therefore we found that the phase model is the strongly interacting system 
and the field theory with non-local interactions on the lattice. 

3.2 Algebraic Bethe Ansatz for the phase model 

In this subsection, we consider the algebraic Bethe Ansatz for the phase model. We follow 
the convention of [UJ. The L-matrix of the phase model at a site n (n — 1, • • • ,N) is 
defined by 



1 W X n 



G End[C 2 (/z)] g> $, 



(3.7) 



where \x G C is a spectral parameter. Here, the L-matrix is a matrix in an auxiliary space 
C 2 . This L-matrix satisfies the Yang-Baxter equation 



R(fjL, v){L{n) <g> L{v)) = {L{v) <g> L(fi))R(fji i v) 



(3i 



with the R-matrix 



R(ji, v) 



(-*- \ 
0-2-1 

fi—V 

0^0 



fj,— V 



-^ 



G End[CV)®C 2 (z/)]. 



(3.9) 



The monodromy matrix is defined by 

T(jm) = L n ^)L n . 1 { / j;)---L 1 (jm) = • (3-10) 

From the Yang-Baxter equation (I3.8p . the monodromy matrix satisfies a following relation: 



R(fi, i/)(T(/x) ® T(z/)) = (T(i/) ® T(/j,))R(n, v). 



(3.11) 
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From this relation, we can derive commutation relations of the monodromy matrix ele- 
ments, A(/j,), B(/jl), C(/i), D{jj): 

B(n)B(v) = B(v)B(n) (3.12) 

C(n)C(v) = C(u)C(fi) (3.13) 

(// - v)A(n)B(v) = -vB{v)A{n) + uB(fi)A(u) (3.14) 

(ft - v)D{n)B{v) = fiB(i/)D(fj) - vB{n)D{v) (3.15) 



v 



C(jm)B(v) = — (A(u)D(jjl) - A(jm)D(u)) (3.16) 

Taking trace with respect to the auxiliary space, the monodromy matrix becomes the 
transfer matrix 

T(»)=tTT(») = A(ri + D(ri. (3.17) 

The vacuum state |0) and its dual vacuum state (0| satisfy 

C(/i)|0) = and (0\B(fj) = 0, (3.18) 

because C(f/,) is the creation operator and B(/i) is the annihilation operator. a(/x) = 1 
and d(fi) = /j n are the eigenvalues of operators A(fx) and D(fi) on the vacuum state |0), 
respectively. 

Eigenstates of the transfer matrix can be constructed by repeated actions of the oper- 
ators B(X) on the vacuum state |0), that is to say, a state fljLi ( 1 ) 1 0} is the eigenstate 
of the transfer matrix 

k k 

r^)f[B{K l m = A(/i,{A})n5(A; 1 )|0) (3.19) 

a=l a=l 

where 

a(m,w) =n_L- +M ''n-i% T , (3.20) 

a=l o=l ^ a 

if the spectrum parameters {Aj}^^...^ satisfy Bethe Ansatz equations 

Af+^-n A," 1 = 1. (3-21) 

6=1 

Finally, we consider a Bethe norm which is a norm of the Bethe vector ^({A}^)) = 
rC=iS(V)|0> and its dual vector (^({A}*)| = (0| flti ^A" 1 ). Here, {A, } j=1 ,., fe satisfy 
the Bethe Ansatz equations. The Bethe norm in the phase model is given by 

mMk)m{\}k)) = tJ^^-mn + kf-\ (3.22) 
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See [9], [16] for the derivation. In the next section, we will clarify connections between the 
U(N)/U(N) gauged WZW model and the phase model. 

4 The U(N) /U(N) gauged WZW model and the phase 
model 

In this section, we clarify connections between the U(N)/U(N) gauged WZW model and 
the phase model. First of all we have to identify parameters of both theories. We identify 
the level k and the rank iV of the gauge group U(N) in the U(N)/U(N) gauged WZW 
model with the total particle number k and the total site number N in the phase model, 
respectively. Under these parameter identifications, we can show that the constraints 
(12.171) coincide with the Bethe Ansatz equations (13.211) . Taking the parameterization of 
the Bethe roots as A a = e 2nl ^ a , the logarithm form of the Bethe Ansatz equations becomes 

{N + k)()) a -J2<l ) b + ^- = m a (4.1) 

6=1 

where m a G Z implies branches of the logarithm. Once we identify the constant field 
(p a in the U(N)/U(N) gauged WZW model with the Bethe roots a in the phase model, 
we found that these equations coincide with the localization configurations (12.171) in the 
U(N)/U(N) gauged WZW model. 

Next, let us consider solutions of the Bethe Ansatz equations. The solutions are ( 12.181) 
because the Bethe Ansatz equations are equal to the localized configurations (12. 17ft . Then 
we can show that piecewise independent solutions of the Bethe Ansatz equations coincide 
with the solutions to be included in the range of < <p a < 1 and to satisfy the condition 
< 0i < 02 < • • ■ < <Pk in the U(N)/U(N) gauged WZW model. Thus, we found that 
this solutions of the Bethe Ansatz equations coincide with {Sol}. The solutions (12.181) 
also imply the completeness of the state in the phase model because the number of the 
solutions is {N + k - 1)!/(JV - l)!Jfe!. 

Since the Bethe norm in the phase model ( I3.22p becomes 

1 r k 2 , Ki(k— 1)0 

me 2m %)m{e 2m %)) = "T 1 * z—Ak + N^N (4.2) 

under taking the parameterization of the Bethe roots as A a = e 2m< ^ a , the partition function 
( I2.23P of the U(N)/U(N) gauged WZW model can be represented as 

« } w(^)= E me 2m %m{e 2 ^} k )) h - 1 . (4.3) 
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Why the partition function of the U(N) /U(N) gauged WZW model can be represented 
by the Bethe norm in the phase model? To understand this, we recall that the partition 
function is represented by using the modular S-matrix ( 12.211) . Thus, we can expect that 
there is a relation between the modular S-matrix in U(N)/U(N) gauged WZW model 
and the Bethe norm in the phase model. Actually, Korff and Stroppel constructed the 
Verlinde algebra in the SU (N) WZW model on the sphere from a viewpoint of the phase 
model and showed that the modular S-matrix in SU(N) WZW model coincides with 
the Bethe norm in the phase model [TT]. So, let us derive the partition function of the 
SU(N)/SU(N) gauged WZW model from the one of the U(N)/U(N) case. There are 
two differences between these partition functions. Firstly, the modular S-matrices in each 
model are related to 



q u(N) _ j N 9 iu(iV) u 4 x 

^oiz -VJ^Tk^oR ■ {lA > 



where 1Z and R denote the primary field in the U(N) and the SU(N) WZW model, 
respectively [15J. Secondary, a range which the summation runs through is different 
because the number of the each primary field is different. Taking account these two 
differences, we find that the partition function of the SU (N) / SU (N) gauged WZW model 
is 

i n^=i(e 27r ^-e 2 ^) Nl ^ 



,<fee{Soi} 



and can be represented by the summation of the Bethe norm with respect to all the 
eigenstates of the transfer matrix in the phase model; 

2S(S,)= £ m{e^h)m{e™%)) h -\ (4.6) 
0i,-,0fce{Soi} 

This shows that the modular S-matrix of the SU(N) WZW model coincides with the 
Bethe norm; 

Sm N) = me 2m *}kM({e 2m nk)). (4.7) 

This is considered as a reason why the partition function of the U(N)/U(N) gauged WZW 
model can be represented by the Bethe norm in the phase model. Therefore, we found 
that the Gauge/Bethe correspondence between the G/G gauged WZW model and the 
phase model is also considered as the gauged WZW model realization of [TT] . 
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Finally, we comment relations between the CS theory and the phase model. The 
partition function of the G/G gauged WZW model coincides with the partition function of 
the CS theory with the gauge group G on S 1 x [6]. We can apply equivariant localization 
methods to the CS theory in a similar way with G/G gauged WZW model. Thus in the 
CS theory with the gauge group U (N) on S 1 xT, h , the localization configurations coincide 
with the Bethe Ansatz equations and the partition function is represented by the Bethe 
norm in the phase model: 



Z^iS 1 x E h ) = 




Further, when the gauge group is SU(N), the partition function of the CS theory is 

ZS* W ( S 1 xEh)= (H{e 2m ' p } k )m{e 2m4 >} k )) h - 1 . (4.9) 

<fr.,-,<frte{Soi} 

We have shown the Gauge/Bethe correspondence between CS theory on S* 1 x S/j and 
the phase model. The equivariant localization for the CS theory on wider class manifolds 
(Seifert manifolds) is derived in [17], [18] and [19], see also [20] for generalization to the 
Chern-Simons-Matter theories. To describe the partition function of the CS theory on 
these manifolds, not only modular S-matrix but also modular T-matrix is needed. It 
would be interesting to consider meaning of T-matrix in the phase model side. 

5 Conclusion 

In this paper, we have studied the relation between the U(N)/U(N) gauged WZW model 
and the algebraic Bethe Ansatz for the phase model. We found that the localization 
configurations ( 12. 17ft coincide with the Bethe Ansatz equations ( 14. ip . once the diagonal 
group elements, the level and the rank of the gauge group U(N) in the U(N)/U(N) 
gauged WZW model are identified with the Bethe roots, the total site number and the 
total particle number in the phase model, respectively. We also showed that the partition 
function of the U(N)/U(N) and the SU(N)/SU(N) gauged WZW model is represented 
as the summation of the Bethe norm with respect to the all eigenstates of the transfer 
matrix in the phase model. This is because the modular S-matrix in the SU(N) WZW 
model coincides with the Bethe norm. This is also considered as the gauged WZW model 
realization involving a generalization to a higher genus case of [IT] . We further found that 
the partition function of the CS theory on S 1 x is also related to norms of Hamiltonian 
eigenstates for the phase model. These relations are summarized in the table 1. 
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the phase model 


the U(N)/U(N) GWZW model/the U(N) CS theory 


Bethe root 


diagonal group element /holonomy along S 1 direction 


Bethe Ansatz equation 


configuration of (J2.17J) 


Total site number 


rank of the gauge group U(N) 


Total particle number 


level 


Bethe norm 


modular S-matrix Sop, 



Table 1: Dictionary in the Gauge/Bethe correspondence 

We comment on some directions for future works. We have considered the case of 
the gauge group U(N) or SU(N) in the G/G gauged WZW model. It is interesting to 
consider a generalization to an arbitrary semi-simple gauge group G. 

Equivariant localization formula itself is not restricted to the topological Yang-Mills 
theory [21] nor the G/G gauged WZW model. One can also couple additional matter fields 
to the topological Yang-Mills theory and the G/G gauged WZW model, and perform the 
localization. Then, these systems will correspond to other integrable systems. For an 
example, the localization configuration for the G/G gauged WZW model coupled to a 
one-form Higgs field leads to the Bethe Ansatz equations for an infinite spin limit for 
higher spin XXZ model |2J. The relation between the partition function of the G/G 
gauged WZW + Higgs system and the norm of Bethe vector for XXZ model will be 
studied in detail [22] , 

Another example is following: Nekrasov and Shatashvili found that a twisted superpo- 
tential in M = (2,2) supersymmetric gauge theory coincides with a Yang- Yang function 
for corresponding integrable system [23], [21]. When we consider the topologically twisted 
supersymmetric gauge theory on a compact Riemann surface with same matter content 
of [23], [24], the localization configuration will lead to the same Bethe Ansatz equation 
which arises from the saddle point equations for the effective twisted superpotential of 
M = (2, 2) theory. For the localization of A/" = (2, 2) supersymmetric gauge theory on S 2 , 
see [25]; [26J. It would be interesting to study localization for such topological twisted 
gauge theories and correspondence with integrable systems. 
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